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ABSTRACT 

Standard two-parameter compressions of the infinite dimensional dark energy model space show 
crippling limitations even with current SN-Ia data. Firstly they cannot cope with rapid evolution - 
our best-fit to the latest SN-Ia data shows late and very rapid evolution to wq = —2.85. However 
all of the standard parametrisations (incorrectly) claim that this best-fit is ruled out at more than 
2<j primarily because they track it well only at very low redshift, z < 0.2. Further they incorrectly 
rule out the observationally compatible region w <^ —1 for z > 1. Secondly the parametrisations give 
wildly different estimates for the redshift of acceleration, which vary from Zacc = 0-14 to Zacc = 0.59. 
Although these failings are largely cured by including higher-order terms (> 3 parameters) this results 
in new degeneracies and opens up large regions of previously ruled-out parameter space. Finally we 
test the parametrisations against a suite of theoretical quintessence models. The widely used linear 
expansion in z is generally the worst, with errors of up to 10% at z = 1 and 20% at z > 2. All of this 
casts serious doubt on the usefulness of the standard two-parameter compressions in the coming era 
of high-precision dark energy cosmology and emphasises the need for decorrelated compressions with 
at least three parameters. 
Subject headings: 



1. INTRODUCTION 

The issue of dark energy dynamics is perhaps the most 
pressing today in cosmology. There are claims both 
for and against dynamics (Bassett et al 2002, Alam 
et al. 2003, Jassal et al. 2004) . But it is a sub- 
ject dogged by gau^e problems ^ (iJonsson et al. 2004k 
IVirev et al. 2004HWan g and Tegmark 2004\ 

Claims for dynamics (or the lack thereof) can be pure 
gauge artifacts, mirages induced by the parametrisations 
that are unrelated to the data. For example, Riess et 
al. (2004) and Jassal et al. (2004) claim that current 
SN-Ia data are inconsistent with rapid evolution of dark 
energy. Such conclusions must always implicitely refer 
to a finite dimensional subspace of the full dark energy 
model space, and broadening the class of models studied 
can (and in this case does) lead to complete reversal of 
such conclusions. Fig. ([2]) provides a explicit counter- 
example. 

The aim of compression is to summarise, as aggres- 
sively as possible, the key features of dark energy prop- 
erties in a few parameters and to facilitate discrimination 
between models. We will see that standard expansions 

^ Consider two dark energy parametrisations. Pi and P2. We are 
interested in estimating the true value of an observable O, such as 
w(z = 0). Under a change in parametrisation the best estimate of 
O will change to O2 = Oi -\- SO. This is a gauge artifact since the 
change in O has nothing to do with the real universe. Even worse, 
since the two parametrisations will usually be rather different we 
have no guarantee that 60 will be small. 



fail to achieve either of these goals, even when they are 
extended to higher order. 

The main result of this work is that compression of the 
dark energy space into low-dimensional subspaces, while 
convenient and easy to work with, can give seriously mis- 
leading results. Since these results are used in the design 
of upcoming surveys it is bad news for cosmology in gen- 
eral. If one does not impose the weak energy condition 
(WEC), w > —1, then the results can border on the 
completely useless. The rest of this article delimits, as 
precisely as possible, the quicksands and danger areas in 
the use of two-parameter compressions. 

As a first sobering example, consider constraints on 
w{z) when we do not impose the WEC. One-parameter 
studies give constraints such as —1.38 < wde < —0.82 
at 2(7 (iMelchiorri et al. 2003f ) suggesting that a model 
with w = —5 at 2; = 2 would be ruled out at more 
than lOcr. Instead a little thought makes it clear that 
if w{z) can vary freely then there is no lower bound 
on w for z > 1 since this merely changes how fast 
the already irrelevant and rapidly diminishing dark en- 
ergy density decreases. If the rapid drop in w oc- 
curs a t z > 1 this leaves essentially no obs ervable 
trace ( Bassett et al. 20021 : ICorasaniti et al. 2003 t). This 
is clearly reflected in the likelihoods in Fig. ([T]) which 
allow for < —100 at z ~ 1. How can we hope to cover 
such possibilities with simple one or two parameter com- 
pressions? 

The dark energy literature overflows with one. 
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two and higher-dimensional compressions of wde{z)^ 
e.g. (Efstathiou 1999; Huterer and Turner 2001; 
Weller and Albrecht 2002; Bassett et al. 2002; 
Corasaniti and Copeland 2002; Linder 2003; Jassal 
et al. 2 004). Comp ressions also exist for Pde{z) 
(|Wang and Freese 20041 : IWetterich 2004 ) while decorre- 
l ated reconstructions of wpEi z) ha ve been p roposed in 
(|Huterer and Starkmann 2003[ ) and (|Hu 200 4'). 

The most basic parametrisation, namely describing the 
dark energy with a constant equation of state w = wq^ 
is well-known to introduce a severe bias (see for instance 
Maor et al 2002, Virey et al 2004) in parameter es- 
timation. Compressions invoking two parameters which 
somewhat alleviate this problem have been introd uced in 
dEfstathiou 1999"; 'Linder 2003"; [JassaLelaL 2004) . How- 
ever, as we will see, these models all struggle to describe 
rapid evolution. This is not surprising. With two pa- 
rameters one may fix at z = and w at high but 
one can do nothing about the time nor the rapidity of 
the transition between the two extremes. Caldwell and 
Doran (2004) circumvented this by considering thirteen 
different one and two-parameter models, some exhibiting 
rapid transitions. 

The use of more than two parameters offers the op- 
portunity to test the dark energy evolution with differ- 
ent data sets and consistently account for the dark en- 
ergy perturbations at all redshifts (Bassett et al. 2002; 
Corasaniti and Copeland 2002; Corasaniti et al. 2004) 
but is, of course, more computationally intensive, and 
care must be taken to accurately capture the dark en- 
ergy dynamics of rapid transitions; see Appendix C of 
(Corasaniti et al. 2004). 

2. THE PARAMETRISATIONS 

For our study we consider two distinct classes of 
compressions. First are standard Taylor expansions of 
wde{z) and second is the Kink, a physically- motivated 
compression. The Taylor expansions are all of the form: 



(1) 



n=0 



where we consider four different choices for the 'expan- 
sion' functions, Xn{z). Namely: 



xo{z) 

Xn{z) 



1 ; Xfi 



0, n > 1 



Xn{z) = {l-ar = 

Xn{z) = {\0g{l ^ Z)) 



l^z 



(constant w) (2) 
(redshift) (3) 

(scale factor) (4) 

(logarithmic) (5) 



To linear order (n < 1) these were first discussed by 
Huterer and Turner (2001) & Weller and Albrecht (2002), 
Chevallier and Polarski (2000) and Linder (2003), and 
Efstathiou (1999) for the redshift, scale-factor and loga- 
rithmic expansion functions respectively. Later we will 
consider their performance at higher order (n > 2). 

The Kink, on the other hand, is not an expansion. It is 
a 4-parameter model which accurately captures the be- 
haviour of quintessence (Bassett et al. 2002, Corasaniti 
and Copeland 2002, Corasaniti et al. 2004). The extra 
parameters allow us to model very rapid transitions in 
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Fig. 1. — Expansion order is more important than 
parametrisation. 1-d marginalised likelihoods for the different 
classes of parametrisations. The narrow curves are the linear like- 
lihoods (constant W£)e^ n = 0, is dotted, first-order expansions, 
n < 1, are dashed). The solid lines show the likelihoods at second- 
order (n < 2) and are much wider due to degeneracies and the 
absence of a lower-bound on w for z > 1. The curves correspond to 
expansions in (a) scale-factor, (b) log and (c) redshift respectively. 
Finally the 1-d likelihood for Wm of the Kink is shown (dot-dashed) 
and exhibits strongly non-Gaussian wings that extend to very large 
values of liUmI- For definition of the order of the expansion see eq. 
P- 

u)de{z)^ a freedom we will need: 
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(6) 



a is the scale factor, wq and Wm are the present and 
matter-dominated values of the dark energy equation of 
state, at is the value of the scale factor at the transition 
from Wm to Wo and A controls the width of the transition. 
Other formulations of the Kink, with relative merits, are 
discussed in Appendix A of Corasaniti et al. (2004). 

There are other parametrisations but these are the 
most widely used today and lessons learned from these 
compressions will apply to many others in the literature. 

3. CONSTRAINTS FROM SN-IA 

In this section we investigate whether the different 
parametrisations above give rise to different best-fits to 
current type la supernova (SN-Ia) data. Recently 6 new 
SN-Ia at redshift z > 1.25 have been discovered using the 
Hubble Space Telescope, providing further evidence for 
a transition from decelerated to accelerated expansion in 
the past(Riess et al. 20 04). 

In or der to be conserva tive we use only the gold sam- 
ple of (|Riess et al. 200^ . containing 157 data points. 
In our analysis we assume a flat Friedmann-Lemaitre- 
Robertson- Walker (FLRW) universe. The assumption of 
flatness is required to achieve reasonable error bars (c.f. 
Kunz and Bassett 2004; Dicus and Repko 2004). Fortu- 
nately this is now a data-driven assumption and particu- 
larly harmless for this study since we are primarily inter- 
ested in testing compressions rather than deriving con- 
straints. We will also assume the prior = 0.27 ±0.04. 
This can be justified from CMB data, as shown in Kunz 
et al. (2003) and Corasaniti et al. (2004), the best-fit 
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values for background FLRW parameters are not affected 
strongly by dark energy dynamics. 
The luminosity distance is given by 

and 

H{Z) = Ho[nm{l + Zf + (1 - nm)f{z)]^^^ (8) 

where f{z) is the solution to the continuity equation 

/(.)=exp(3^^i±i^d.'). (9) 

Our analysis methods are described in detail in Corasan- 
iti et al (2004). We use a Markov-Chain Monte Carlo 
code to find the constraints on the dark energy param- 
eters for each parametrisation. As usual we marginalise 
analytically over the normalisation ofdL which takes care 
of the Hubble constant as well, leaving as the only re- 
maining parameter apart from those describing the equa- 
tion of state. 

Figure [T] shows the marginalised one-dimensional like- 
lihoods for the parameters of the dark energy compres- 
sions. The main point of that figure is that the various 
parametrisations have similar likelihoods at the same or- 
der^ but that the likelihoods at different orders are com- 
pletely different. Here by order we mean the maximum 
value of n in eq. ([1]). Hence "linear" or "first-order" 
(n < 1) refers to the standard two-parameter expansions 
with only wq^Wi non-zero, while 2nd-order corresponds 
to n < 2 and has non-zero 

The best-fit values and la errors for the first-order 
redshift, scale factor and logarithmic parametrisations 
are given in table [U these are consistent with those found 
in (Feng et al. 2004; Gong 2004). 

TABLE 1 

Best fit values and la confidence intervals at linear 

ORDER {n < 1). 



wo Wl 



Redshift 


-i-3o±g:i 


i.57±i:tf 


Scale-factor 


-i-48±C:iI 


3.ii±i:?i 


Log 


-i-39±g:i? 





We infer the 'maximised' limits on the redshift depen- 
dence of the equation of state by computing for a given 
parametrisation the highest and lowest wde{z) for the 
models in the chains with < Xmin + ^- Models with 
w{z) outside those limits are therefore expected to be 
"bad" fits to the SN-Ia data. We plot the result in Fig- 
ure [21 For the two-parameter models, the limits are very 
similar to those obtained by marginalising over all other 
parameters, as is expected since they are nearly Gaus- 
sian. We have also checked that they coincide with limits 
from Gaussian error propagation. 

We now discuss the constraints derived from the Kink 
formula, Eq. (j6|). The best fit to the data has a = 
172.6 and is characterized by: wq = —2.85, Wm = —0.41, 




z 



Fig. 2. — Parametrisations struggle with rapid evolu- 
tion. Maximised limits on wde(z) for the redshift (red dashed 
hne), scale-factor (green dash-dotted), logarithmic (blue dotted) 
and Kink (solid black line) parametrisations. The best-fit kink 
solution passes well outside the limits of all the parametrisations 
(except for the kink) both at z ~ and at z ^ 0.2, showing their 
inability to capture rapid dynamics which leads to their incorrectly 
ruling it out. 

at = 0.94 and log(A) = —1.52, corresponding to a 
rapidly varying equation of state with a transition from 
Wm to Wq at Zt = 0.1. This best-fit is shown in fig- 
ure [2] and clearly exits the 2a limits from the two- 
parameter compressions, first from below at 2: ~ and 
then from above, at z ~ 0.2. This graphically illus- 
trates the limitations of the standard parametrisations 
and shows how they artificially rule out models which 
sho uld give the strongest s ignals for dark energy dynam- 
ics (jCorasaniti et al 20031 ). 

Although the Kink formula is by construction every- 
where under control, the dependence on its parameters 
is highly non-linear so the peaks of the marginalised 
likelihood distributions do not match the values of the 
best fit model. In fact the peaks of the 1-dimensional 
marginalised likelihoods are shifted with the respect 
to the best fit values after marginalising. Hence the 
marginalised likelihoods do not coincide with the maxi- 
mized ones, as they would if the likelihood distributions 
were Gaussian. 

It is not just one good fit to the data which violates the 
2a limits of all the two-parameter compressions either. 
For instance, the model with wq = —1.46, Wm = 0.16, 
at = 0.88 and log (A) = -0.7 has = 173.9, while 
the model with wq = —1.11, Wm = 6.13, at = 0.40 and 
log (A) = -.98, has = 175.9. Both are excellent fits to 
the data but are supposedly ruled out by the n < 1, linear 
redshift, scale-factor and logarithmic parametrisations of 
equations ([3]|5]). 

The conclusion that rapid evolution of dark energy is 
ruled out by current data is therefore a 'gauge' artifact. 
We have shown that rapid variations of the dark energy 
equation of state are perfectly consistent with, and in 
fact provide better fits to the Gold sample than do mod- 
els without rapid transitions. This conclusion remains 
even after including CMB and large scale structure data 
(Corasaniti et al 2004). 
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The pathological behaviour of ruling out models which 
are very good fits to the data can be rectified by the 
inclusion of higher order terms, n > 2. Indeed, since the 
data allow wi to be large in all cases, higher order terms 
in the redshift, scale-factor and Log expansions cannot 
be neglected. Therefore we have extended our analysis 
in order to include second order corrections (n < 2) to 
Eq. da])-©. 

Comparing the yellow likelihoods in Figure [1] with the 
red ones we see that by the allowed values of wq are 
shifted significantly towards more negative values, con- 
sistent with, but broader than the kink confidence inter- 
val. Secondly huge values of wi ~ 50 and W2 ~ —100 are 
consistent with the data. This suggests that strong dark 
energy dynamics is not ruled out and that higher order 
terms must be taken into account. As mentioned in the 
introduction, this comes from the fact that Wi and W2 
are strongly degenerate in all cases and that there is no 
lower bound on il' at z > 1, illustrating the huge effect 
of imposing the weak energy condition w > —1. 

In this case we have both a lower bound (from the weak 
energy condition) and an upper bound (from the data) on 
w. As the second order parametrisations are described by 
a parabola in their respective expansion variables, they 
end up being strongly constrained. 

On the other hand, we have also considered much 
higher order (n < 6), and found that severe internal 
degeneracies lead to finely balanced coefficients. Thus 
a Taylor expansion in z becomes very unstable at high 
redshift, and even expansions in 1— a can hardly be called 
"under control" anywhere since the expansion coefficients 
become of order 10^. 

3.1. Model selection with Information Criteria and the 
Bayesian Evidence 

In table [2] we report the values of the foi" the best fit 
models of each parametrisation. As we can see the four- 
parameter Kink formula has the lowest followed by 
the scale-factor, the logarithmic and redshift parametri- 
sation. Accounting for second order terms (n < 2) in 
the expansions provide fits even better than the best 
Kink parametrisation which is non-trivial since they have 
one less parameter, although now the allowed parameter 
space volume is huge. 

At this point we can ask how many parameters are 
actually necessary to describe the dark energy with cur- 
rent SN-Ia data? Following Liddle (2004), we compute 
the Akaike information criterion (Akaike 1974) 

AIC = -2\nC^2k (10) 

and Bayesian information criterion (Schwarz 1978) 

BIC = -2\nC^k\nN, (11) 

where C is the maximum likelihood, k is the number of 
model parameters and N is the number of data points. 
We also compute the Bayesian evidence E (Sivia 1996; 
Saini, Weller and Bridle 2004), both with thermody- 
namic integration and nested sampling (Skilling 2004). 
It is worth noting that fully degenerate parameters do 
not contribute to the Evidence, so that specifically the 
Kink model is less disfavored than the number of param- 
eters suggests naively. For the same reason we find that 
E grows very slowly when going to even higher order 
in the expansion-type parametrisations, although these 



TABLE 2 

Bayesian Evidence, BIC, AIC and values of best fit for 

DIFFERENT PARAMETRISATIONS. k IS THE TOTAL NUMBER OF 
FITTING PARAMETERS WHICH INCLUDE THE DARK ENERGY DENSITY 



model 


k 


x' 


BIC 


AIC 


-HE) 


ACDM 


1 


177.6 


182.7 


179.6 


93 


w=const 


2 


177.6 


187.7 


181.6 


96 


linear 


3 


174.5 


189.7 


180.5 


99 


logarithmic 


3 


174.2 


189.4 


180.2 


98 


scale factor 


3 


174.0 


189.2 


180.0 


98 


quadratic 


4 


172.1 


192.3 


180.1 


100 


logarithmic II 


4 


172.2 


192.4 


180.2 


100 


scale factor II 


4 


172.3 


192.5 


180.3 


99 


Kink 


5 


172.6 


197.9 


182.6 


96 



TABLE 3 

Best fit values and la confidence intervals on Zacc at 

LINEAR ORDER (n < 1). 



Model 


Zacc 


Model 


Zacc 


ACDM 


0.66±0:}} 


w = const. 


U.CiD^Q Q4 


Linear in z 




Log 




Scale-Factor 


o.59±g:l} 


Kink 


n 4c_l0.53 



cases are already disfavored by Bayesian statistics. The 
preferred parametrisation is ACDM - it is indeed remark- 
able that a model with a single free parameter fits the 
data so well. 

4. WHEN DID ACCELERATION BEGIN? 

One of the key characteristics of dynamical dark en- 
ergy is that the redshift at which the universe begins ac- 
celerating, Zacc^ is characteristically different from that 
in ACDM with the same Qde today. This is manifest 
i n the CMB as a modified integrated Sachs - Wolfe effect 
(Bassett et al. 2002J; ICorasa niti et al. 2003^ which is de- 
generate with reionisation (Corasaniti et al. 2004). The 
SN-Ia offer the possibility to break this degeneracy and 
therefore it is crucial to use a parametrisation that can 
accurately estimate Zacc without bias. Using a simple 
linear expansion of the deceleration parameter Riess 
et al. (2004) estimated Zacc — 0.46 ±0.13. Here we com- 
pare the predictions of the various parametrisations for 
Zacc The Id likelihoods for Zacc are shown in Figure [3] 
and in table [3] we report the confidence intervals. 

We have two main points. First, all of the parametri- 
sations predict different best-fit values and la error bars 
for Zacc ranging from Zacc = 0.14 for the redshift ex- 
pansion to Zacc = 0-59 for the scale-factor expansion (see 
also Dicus and Repko 2004). The logarithmic, constant 
and Kink parametrisations all have similar best-fits, but 
the first two have overly narrow error bars relative to the 
Kink predictions. The largest error bars correspond to 
the scale- factor expansion, eq. (|4]). 

Given the importance of accurately estimating Zacc this 
variance forcefully argues for the need to go beyond two- 
dimensional parametrisations in handling future, high- 
quality data. 

As a second point it is interesting that the best-fits 
for Zacc for all the parametrisations are lower than in 
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the ACDM model. This may be novel evidence for dark 
energy dynamics. We leave this issue for future work. 




Fig. 3. — When did acceleration begin? The marginalised Id 
likelihood for Zacc as estimated from the same SN-Ia data for each 
of the various dark energy parametrisations. The wide variance 
casts doubt on their usefulness. All of them peak at a lower Zacc 
than in the standard ACDM, perhaps giving novel evidence for 
dark energy dynamics. Only the scale-factor expansion, eq. (0) 
shares significant overlap with ACDM, but that is expected since 
the transition is typically slow in that parametrisation. 



5. CONCLUSIONS 

This Letter shows the limitations of standard one 
and two-parameter compressions of the infinite dimen- 
sional space of dark energy models. We have high- 
lighted the dangers in using constraints derived using 
these parametrisations, particularly regarding the possi- 
bility of rapid evolution in the dark energy, which none 
of the standard compressions can follow, and in defin- 
ing allowed regions of parameter space which are com- 
pletely wrong in the case where the weak energy condi- 
tion {w > —1) is not imposed. 

Rapid evolution provides a superlative fit to cur- 
rent SN-Ia data (as measured by x^)? despite claims 
to the contrary in the literature which were based on 
two-parameter compressions. Indeed, all of the two- 
parameter expansions we studied wrongly rule out such 
rapid evolution at 2a or more. This is extremely damn- 
ing evidence, especially since these compressions are typ- 
ically used in the planning for the next-generation ex- 



periments which will provide data of significantly higher 
quality. In addition, the standard parametrisations also 
miss the fact that w has no lower-bound at z > 1 if the 
weak energy condition is not imposed, artifically cutting- 
out vast swaths of parameter space due to their innate 
limitations. 

Further problems occur in estimating the redshift at 
which the universe began accelerating, Zacc- There is a 
nearly 300% variation in the best-fit for Zacc^ depend- 
ing on parametrisation. Interestingly, all the tested 
parametrisations gave best-fits for Zacc below that of 
ACDM, providing unusual cross-parametrisation evi- 
dence for dark energy dynamics. Nevertheless, use of the 
Bayesian information criteria for model selection prefers 
the cosmological constant over the other models which 
remains the model to beat. 

The severe inadequacy of the standard two-parameter 
expansions lead us to consider higher-order terms (n > 2) 
with one of more extra parameters, e.g. While this 
brings the rapid evolution models within the allowed re- 
gion of parameter space it leads to severe degeneracies 
(see figure [T]) which may make the parametrisations im- 
potent for constraining the space of theoretical dark en- 
ergy models, particularly when w < —1. 

For completeness we have also compared four standard 
parametrisations against a test-bed of quintessence mod- 
els (Appendix) and found that while the kink and scale- 
factor expansion are excellent, the expansions in z and 
log(l + z) can lead to large errors for z > 1. 

We conclude that confidence intervals inferred from 
standard two-parameter expansions often do not deserve 
that name and are typically untrustworthy, even with 
current data. The wealth and quality of dark energy 
data we will aquire over the next decade will demand 
significantly better performance. 
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APPENDIX 

In this appendix we address the problem of how accurately the various parametrisations listed earlier do in capturing 
the dynamics of standard quintessence models. We proceed by considering a test-bed of four popular quintessence 
models, which provide a wide range of different behaviours for w{z), with smooth, slow and rapid transitions. Our 
test-bed suite of models is: the Albrecht-Skordis model (Albrecht and Skordis 2001), 

V{Q) = M^e-"^[(g - Af + (1) 

with a = 5, A = 54.5, B = 0.01 and s = 5; the Two Exponential potential (Barreiro, Copeland and Nunes 2000), 

V(Q) = M'^(e-"^ + e-^^), (2) 

with slopes a = —9 and /3 = 2; the inverse power-law potential (Ratra and Peebles 1988), 

ViQ) = M^+^/Q", (3) 

with a = 0.6; the Supergravity (SUGRA) inspired model (Brax and Martin 1999), 

V{Q) = M4e-'5V2/ga^ (4) 
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max max 



Fig. 4. — Plot showing the mean quadratic error of the best-fit, eq. (|5]), to the different quintessence models as function of Zmax for each 
parametrisation for n < 1. At z ~ 1 the expansion to first order in z ('redshift') and log(l + z) ('logarithmic') typically show errors around 
10% while the kink and scale-factor parametrisations have errors typically at the 1% level at z ~ 1. 

with a = 6. For the kink and each two-parameter expansion (n < 1) we compute the best-fit (using an MCMC search) 
to the different quintessence models in the redshift range < z < Zmax and plot the average residual quadratic errors 
as a function of Zraax^ as shown in Figured! where 

1 

= aT E(^Q(^^) - (^0)', (5) 

^ i=l 

with Nz the total number of equally spaced bins in the range [0, Zmax]^ and we have chosen Nz = 1000. Comparing the 
residuals indicates how well each formula reproduces a given quintessence model as function of the redshift interval. 
For Zmax ^0.3, the Kink and scale-factor formulae provide best-fits which are at least one order of magnitude better 
than the redshift and logarithmic expansions, except perhaps in the Two- Exponential case. 

Only for Zmax < 0.3 are the residuals of the redshift and logarithmic formulas of the same order as those of the other 
parametrisations. This is easy to understand from Fig. ([2]). The best-fit solution is almost perfectly linear for z < 0.1 
but then suddenly deviates. A linear fit will be excellent until z ~ 0.1, after which it rapidly becomes bad. 

We do not consider the case of constant w which can at best fit the average value of w{z) and whose errors are 
expected to be the worst. 
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